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Abstract 

Neutrinoless double beta decay (Oi//3/3) induced by superparticle ex¬ 
change is investigated. Such a supersymmetric (SUSY) mechanism of 
Qvf3f3 decay arises within SUSY theories with R-parity non-conservation 
{lf.j). We consider the minimal supersymmetric standard model (MSSM) 
with explicit .^pterms in the superpotential (.^pMSSM). The decay rate 
for the SUSY mechanism of Dvjfl decay is calculated. Numerical val¬ 
ues for nuclear matrix elements for the experimentally most interesting 
isotopes are calculated within pn-QRPA. Constraints on the .^pMSSM 
parameter space are extracted from current experimental half-life lim¬ 
its. The most stringent limits are derived from data on ^®Ge. R is 
shown that these constraints are more stringent than those from other 
low-energy processes and are competitive to or even more stringent than 
constraints expected from accelerator searches. 


1 Introduction 

In the standard model (SM), since B — L conservation is exact, neutrinoless 
double beta decay (0i//3/3), which violates lepton number by two units, is for¬ 
bidden. On the other hand B — L and L violation is expected in theories 
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beyond the SM. That is why 0i//3/3 decay has long been recognized as a sen¬ 
sitive tool to put theories beyond the SM to the test (for reviews see i. i. 
[^]). A variety of mechanisms which may cause 0z//3/3 have been studied in the 
past. The simplest and the most well-known possibility is via the exchange of a 
Majorana neutrino between the decaying neutrons or due to {B — L)-violating 
right-handed currents. Oz//?/? decay has not yet been seen, but limits on various 
model parameters can be deduced (see p, and references therein) from its 
non-observation. 

Recently impressive progress has been achieved in the experimental inves¬ 
tigation of double beta decay, both in the 2z//3/3 and the Oz//5/3 decay mode 
[i]” fhe near future essential advance in this direction is expected. 

Experimental lower bounds on Oz//3/3-decay half-lives are often represented in 
terms of an upper limit on the Majorana neutrino mass {rriy). At least one 
experiment currently in operation will reach a hnal sensitivity of about {rriu) 
= (0.1 — 0.2) eV and has already pushed the existing limit below 1 eV 0. 

In view of the rising experimental sensitivity it is of great interest to pursue 
a more comprehensive theoretical study of the possible mechanisms of Oz//?/? 
decay. 

In this work we investigate contributions to Oz//3/3 decay within supersym¬ 
metric (SUSY) theories with explicit R-parity breaking. R-parity (Rp) is a 
discrete, multiplicative symmetry dehned as Rp = where 

S', B and L are the spin, the baryon and the lepton quantum number. The 
SM helds, including additional Higgs boson helds appearing in the extended 
gauge models, have Rp = +1 while their superpartners have Rp = —1. This 
symmetry has been imposed on the minimal supersymmetric standard model 
(MSSM) (for a review see ||^) to ensure baryon number [B) and lepton num¬ 
ber (L) conservation. However, neither gauge invariance nor supersymmetry 
require Rp conservation. The question whether or not Rp is a good symmetry 
of the supersymmetric theory is a dynamical problem which might be related 
to more fundamental physics at the Planck scale. In general Rp can be either 
broken explicitly vacuum expectation value of the scalar superpartner of the 

ra. 


Rp-odd isosinglet lepton held 


Supersymmetric models with Rp non-conservation 


have been exten¬ 


sively discussed in the literature not only because of their great theoretical 
interest, but also because they have interesting phenomenological and cosmo¬ 
logical implications. 

Existing constraints on .^pSUSY theories are either direct from collider ex¬ 
periments or indirect from low-energy processes [|T^, matter stability 
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np, UlSf and cosmology We will discuss the former two constraints 

later on in more detail in comparison with the bounds from Oz//3/3 decay. 

Recently .^pSUSY models have been analysed in connection with current 
and forthcoming collider experiments. The .^pSUSY gives rise to spectacular 
signatures of events in collider detectors which would yield a very clean signal 
for supersymmetry. Consequently, expected sensitivities of experiments at 
HERA |2^, the TEVATRON |^|, LEP 200 and the LHC have been 


analysed recently. In the present paper we pay special attention to comparing 
the capability of Oz//3/3 decay experiments with the collider experiments in 
establishing better constraints on .^pSUSY models. 

The rest of the paper is organized as follows. In the next section, we specify 
the minimal supersymmetric standard model with Rp non-conservation. Sec¬ 
tion 3 discusses the basic diagrams inducing decay in the MSSM. 

The effective low-energy Lagrangian as well as the different lepton number vi¬ 
olating parameters are dehned. Section 4 outlines the procedure for obtaining 
nucleon matrix elements from the quark currents in the non - relativistic im¬ 
pulse approximation. Section 5 then deals with the numerical calculation of 
the relevant nuclear structure matrix elements. We briefly summarize the main 
features of the nuclear structure model, before presenting numerical results for 
those isotopes which are currently the experimentally most promising. Special 
attention is paid to a discussion of the theoretical uncertainties of the nuclear 
structure calculation. In section 6, on the basis of the current experimental 
limit on the half-life of ™Ge [||, we analyse constraints on the supersymmet¬ 
ric parameter space imposed by the non-observation of QuPP decay. We have 
found that these limits are more stringent than those from other low-energy 
processes and also more stringent than those expected from experiments with 
the ZEUS detector at HERA We then close with a short summary and 
outlook. 


2 Minimal Supersymmetric Standard Model 
with i?-parity Non-conservation 

In the following we will use the MSSM extended by inclusion of the explicit 
R-parity non-conserving terms {^p) into the superpotential. This model has 
the MSSM held content and is completely specihed by the standard SU{3) x 
SU{2) X U{1) gauge couplings as well as by the low-energy superpotential and 
’’soft” SUSY breaking terms [Q. The most general gauge invariant form of the 
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superpotential is 


( 1 ) 


W = W.+ . 


The Rp conserving part has the standard MSSM form 

Wn = hLHiLE + hoHiQD - huH^QU - /xi/iifa- 


( 2 ) 


We use notations L, Q for lepton and quark doublet superfields and E, f7, D 
for lepton and up, down quark singlet superfields; Hi and H 2 are the Higgs 
doublet superfields with a weak hypercharge Y = —1, +1, respectively. Sum¬ 
mation over generations is implied. For simplicity generation indices of fields 
and Yukawa coupling constants Hl, hu, ho are suppressed. The mass-mixing 
parameter /r is a free parameter describing mixing between the Higgs bosons 
H 1 -H 2 as well as between higgsinos H 1 -H 2 . 

TO 


The Rp violating part of the superpotential (|^) can be written as 


Wu 


^ijkLiLjEj^ -|- jEj« -|- 


(3) 


where indices i,j, k denote generations, and the fields have been defined so that 
the bilinear lepton number violating operators LiH 2 have been rotated 


away. The coupling constants A (A") are antisymmetric in the first (last) two 
indices. The first two terms lead to lepton number violation, while the last 
one violates baryon number conservation. 

Proton stability forbids the simultaneous presence of lepton and baryon 
number violating terms in the superpotential 0 (unless the couplings are 
very small). Therefore, only A, A' or X" type interactions can be present. 
There may exist an underlying discrete symmetry in the theory which allows 
either the first or the second set of couplings |0, An example of such a 
symmetry, which forbids baryon number violating couplings but allows lepton 
violating ones, is given by the transformation rules 


(Q, U,d)-^- (Q, U, d) , (l, E, Hi, 2 ) + (t, E, Hi, 2 ) . 


(4) 


This discrete symmetry can be justified on a more fundamental level of Planck 
scale physics. It has been shown to be compatible with the ordinary SU{5) 
[pud and ’’flipped” SU{5) x 17(1) grand unification (GUT) scenarios as 
well as with phenomenologically viable superstring theories TOl . 

Neutrinoless double beta decay, which is the main subject of the present 
paper, requires lepton number violating interactions. Therefore we bind our¬ 
selves to the ^p MSSM with lepton number violation (A 7 ^ 0, A' 7 ^ 0) and 
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baryon number conservation (A" = 0). The Lagrangian of this model pos¬ 
sesses the discrete symmetry eq. (H). Apparently, Of(3(3 can probe only the 
hrst generation lepton number violating coupling because only the first 
generation fermions u, d, e are involved in this process. 

In addition to proton decay constraints on .^pcouplings there are also con¬ 
straints which follow from cosmological arguments, requiring that the baryon 
asymmetry generated at the GUT scale is not washed out by B — L violat¬ 
ing interactions present in eq. (j^). These cosmological constraints have been 
thought to affect all .^pcouplings A, A', A" -C 10“^, making these models phe¬ 
nomenologically not interesting. These arguments, however, were proved to be 
strongly model dependent m , bounds can be evaded in perfectly reasonable 
scenarios of matter genesis m- 

The effect of ’’soft” supersymmetry breaking can be parametrized at the 
Fermi scale as a part of the scalar potential: 


Vsoft= E rnl\^i\'^ + hLALHiLE + hDAnHiQD-huAuH2QU- (5) 

i=scalars 

—IJ.BH 1 H 2 + h.c. 


and a ’’soft” gaugino mass term 


^GM — -- 


2 L 


M^BB + 


h.c. 


( 6 ) 


As usual, M 3 2,1 are the masses of the SU{3) x SU{2) x U{1) gauginos g, W, B 
and rrii are the masses of scalar helds. Al, Aj:,, Ajj and B are trilinear and 
bilinear ’’soft” supersymmetry breaking parameters. All these quantities are 
free SUSY model parameters which due to the renormalization effect depend 
on the energy scale A. 

Considering a GUT scenario within the MSSM one can claim the following 
unihcation conditions at the GUT scale A ~ Mx'- 


Au{Mx) — Ad{Mx) — Al{Mx) — Aq, (7) 

mL{Mx) = rriEiMx) = mqiMx) = mu{Mx) = moiMx) = mo, (8) 

Mi(Mx) = M2{Mx) = M^{Mx) = mi/2, (9) 

9i{Mx) = = 93i^x) = 9 guti (10) 


where 9 ^, 92,91 are the SU{3) x SU{2) x f/(l) gauge coupling constants equal 
to 9 gut fh® unification scale Mx- 
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At the Fermi scale A ~ Mw these parameters can be evaluated on the basis 
of the MSSM renormalization group equations (RGE) 
that the ^pYukawa coupling constants A, A', A" are small enough to be neglected 
in these equations. Equation @ implies at A ~ Mw 

5 

Ml = - tan^ 6w ■ M2, M 2 — 0.3mg. (11) 

o 

Here rrig = M3 is the gluino mass. 

Now the model is completely specihed and we can deduce the interaction 
terms of the .^pMSSM - Lagrangian relevant for neutrinoless double beta decay. 

Write down these interaction terms explicitly. Note that in the following 
we use for fermion helds the 4-component Dirac bispinor notation. 

The lepton number violating part of the Lagrangian can be obtained di¬ 
rectly from the superpotential (^. It has the form 


. We assume 


(i )*»» 


The Lagrangian terms corresponding to gluino Cg and neutralino inter¬ 
actions with fermions tjj = {u,d,e}, q = {u,d} and their superpartners = 
{u, d, e}, q = {u, d} are ^ 

A(“) 

- QrMr) + h.C; (13) 




A) 


111 


{ul d 


R) 


'R 


'R 


dR + (cl ^l) dji 


+ (ur di) dR ■ 


h.C. 


-V} 


= 1 / 25(2 {eLiWi>LXi^L + + h.c. (14) 

i=l 

Here A*-®^ are 3x3 Gell-Mann matrices (a = 1,...,8). Neutralino coupling 
constants are dehned as [§] 

eLi(V') = -T3{^|J)^^i2 + tanewiT3{^/J) - Q{ilj))Afii, (15) 

(^Riii>) = Q{'i/j)tanewMi- (16) 

Here Qi'ip) and T 3 {'ip) are the electric charge and weak isospin of the held 
i). 


6 






ij are elements of the orthogonal neutralino mixing matrix 


Coefficients 

which diagonalizes the neutralino mass matrix. In the Ijip MSSM the neu 
tralino mass matrix is identical to the MSSM one 
{B, H 2 ) has the form: 


and in the basis of fields 


/ 


Mv = 


Ml 
0 

V MzSy^Sfj 


0 

M 2 

—MzCyYSj^ 


-MzSyyCp 

Mzc-^^CjS 

0 

-/i 


MzSy^Sf} \ 
— MzCy^rSi3 
-/i 

0 y 


(17) 


where = cos9wi s^r = sin^vi/, = tan6'iy, = sin/?, C/? = cos/?. The 

0 


Here < HI > and < Hi > 


angle /? is dehned as tan/? =< H^ > / < H^ > 
are vacuum expectation values of the neutral components H^ and H^ of the 
Higgs doublet helds with weak hypercharges Y{Hl) = +1 and Y{Hl) = —1, 
respectively. The mass parameters Mi, M2 are related to the gluino mass rrig 
according to eq. (Q)- 

By diagonalizing the mass matrix ( p/Tf ) one can obtain four neutralinos Xi 
with masses m^. and the held content 

Xi 


Xi — ■X'iiB + J\fi2W^ + Afi^Hl + AfiiHl- 


(18) 


Recall again that we use notations W^, B for neutral SU{2)l x U{1) gauginos 
and Hi, H^ for higgsinos which are the superpartners of the two neutral Higgs 
boson helds H^ and H^. 

We apply a diagonalization by means of a real orthogonal matrix A/”. There¬ 
fore the coefficients A/ij are real and masses m^. are either positive or negative. 
The sign of the mass coincides with the CP-parity of the corresponding neu¬ 
tralino mass eigenstate Xi- If necessary, a negative mass can be always made 


positive by a redehnition |35] of the neutralino held Xi- H leads to a redehnition 


of the relevant mixing coefficients A/),- i ■ Mi 


‘ly 


The lightest neutralino is commonly assumed to be the lightest supersym¬ 
metric particle (LSP). That is true in almost all phenomenologically viable 
SUSY models with Rp conservation. For .^pSUSY models this is a very non¬ 
trivial assumption, since the cosmological constraints 
to be colour and electrically neutral, no longer apply ^0 
could be any superparticle in ^p SUSY models. However, the RGE analy¬ 
sis in minimal supergravity models suggests that the LSP is a neutralino if 
couplings are reasonably small the lightest neutralino. If Rp is conserved 
the LSP is a stable particle. Otherwise it decays into ordinary matter. 


3g] requiring the LSP 
A priori, the LSP 
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Squarks ul,r, dL,R and selection eL,R in eqs. ([T3|), ([T^) are with a good 
precision mass eigenstates. Possible fi — /rj-mixing for the hrst generation of 
sqnarks and sleptons are negligible dne to the smallness of the relevant Ynkawa 
conplings. In this case the MSSM mass formnlas can be written as 

H 


33 



= ml + 0.07ml+ -cos2pMl{2sm‘^ 9w — 

^ 2 

(19) 


= ml + 0.02m? — cos2/9M| sin^ 9w, 

(20) 


1 4 

= mg-I-0.83m?-I-- cos2/dM|(l — - sin^ ^w), 
^2 3 

(21) 


1 2 

= mg-I-0.83m? — - cos 2/3M|(1 — - sin^ 6 *vk), 
^2 3 

(22) 


2 

= mg -I- 0.77m? -|- cos2/3M|- sin^ 9w, 

^ 3 

(23) 

m} 

= mg-I-0.76m? — cos2/9M|-sin^6*v(/. 

(24) 


Here, rrig is the glnino mass and rriQ is the common sfermion mass at the 
nnihcation scale (see (^)). From the eqs. (0)-(El) one can estimate in the 
region m^j^^^SOOGeV that ~ approximate relation 

will be helpfnl for nnderstanding the resnlts of our numerical analysis in Sect. 

6 . 

Having specihed the Lagrangian interaction terms ([T^) - (|^ and the mass 
eigenstates g, Xi, Q, e involved in the interactions we can construct diagrams 
describing the .^pMSSM contribution to the neutrinoless double beta decay. 

In principle, in the MSSM a small neutrino mass of Majorana type 
may arise radiatively [^, [^. On the other hand in some GUT scenarios 
based on large gauge groups like 50(10) there might be also heavy Majorana 
fermions having non-negligible SU{2)l components. It is often identihed with 
the heavy Majorana neutrino N. In the presence of either the light Majorana 
neutrinos or the heavy ones (or both) Oz//3/3 decay can be induced by the con¬ 
ventional neutrino mass mechanism. In the present paper we concentrate on 
the .^pSUSY mechanism of Oz//3/3 decay. The effect of additional contributions 
from the neutrino mass mechanism is investigated in the case of heavy Majo¬ 
rana neutrino exchange. (Inclusion of the light Majorana neutrino contribution 
will not change the results of our analysis.) 








3 Low-energy ALg = 2 - Effective Lagrangian 


The basic diagram corresponding to the neutrinoless double beta decay at the 
nucleon level is presented in fig. 1(a). Two neutrons from the initial nucleus Ai 
after interaction transform into two protons of the final nucleus Af emitting 
two electrons. Apparently this process breaks the electron lepton number by 
two units, ALg = 2. At the quark level it can be induced by the subprocess 
with two initial d-quarks and two final u-quarks accompanied by two electrons 
as shown in Fig 1(b). 

The conventional mass mechanism with Majorana neutrino exchange is 
presented in fig. 2(a). In the following we consider the .^pMSSM contribu¬ 
tion to Ozz/d/d decay. Starting from the fundamental interactions (B)-(Il) we 
have found the complete set of diagrams presented in fig. 2 (b, c) which 
contribute to this subprocess. 

The supersymmetric mechanism of 0z//3/3 decay was first proposed by Mo- 
hapatra and later studied in more details by Vergados 

0 ^ 


In these papers 

|40[| only three diagrams similar to those in fig.2(b) were considered. In¬ 


stead of neutralinos Xi, which are actual mass eigenstates in the MSSM, the 
consideration of refs. |^, used Z-ino (Z) and photino ( 7 ) fields in in¬ 
termediate states. Z and 7 can be mass eigenstates only at special values 
of parameters of the neutralino mass matrix. In general, these fields are not 
mass eigenstates. Furthermore, it was recently realized that a cosmologically 

0 


viable lightest neutralino is very likely B-ino dominant 


Using in the in¬ 


termediate states such fields which are not mass eigenstates leads to neglecting 
diagrams with mixed intermediate states when, for instance, Z turns to 7 due 
to the mixing proportional to the relevant entry of the neutralino mass matrix. 
The effect of mixing is taken into account completely in the set of diagrams 
displayed in fig. 2 (b,c) with all neutralino mass eigenstates Xi involved. 

In the case of 0i//3/3 decay when momenta of external particles are much 
smaller than intermediate particle masses one can treat interactions in fig. 
1(b) and fig. 2 as point-like. A suitable formalism in this case is the effective 
Lagrangian approach. 


Define the effective Lagrangian £ 


AZ/e — 2 1 
eff ' 


x) as 


</|A-l|*> = 


= i J d^x < f\£e// ‘^{x)\i > + 

high orders of perturbation theory. 


(25) 


Thus, £ef^f “^{x) corresponds to the lowest order operator structure having 
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non-vanishing matrix elements of the form 


^ I /* AZ/e = 2 i 

< uuee\£^fj ( 


x)\dd > 7 ^ 0. 


(26) 


It is now straightforward to hnd the operators in the effective Lagrangian 
which correspond to the diagrams in £g. 2(b,c). The result is 


pALe =2 

^eff 


4 

(x) = SvrasA'niXl 


1- i my. 

2 — J. ACi L 


^^^(u2dRa)(u^djif3)(eLe%)+ (27) 


+ 


'-'Li 


» 


m. 


«2wR/3)(eLC^Ra)(eL4) + 




eL 

(d)^ 


Ke]i)(uLe%)(d^LadRf3) + 


UL 


4 

mA 

dR 


eLi{u)eRi{d) eLi{u)eLi{e) eLi{e)eRi{d) 

“T O 0 “T 


{uldR){uL^e%){eLdRa 


2 2 

dR 




+ Aiii 


Sttck^ '■ 


m-n 


2 2 


2 2 
mi m r 

ei dfl 


«Wf7)(eL4)(eLO + 


+ 




4 

mj 

dR 




2 2 
mj mj. 


{"^I^r) i^id^) 


Gauge coupling constants Q (2 = 5'|/(47r) and = 5f|/(47r) are running coupling 
constants which should be estimated at the proper energy scale A. One can 
see from the diagrams in £g.2 that the typical scale at which the g — q — q and 
X~q~q interactions occur is of the order of the gluino or neutralino mass. In 
the mass region which will be analysed numerically in the subsequent sections 
we may take approximately these couplings at the Z-boson pole and use their 
values from ref. |^ : 


as{Mz) = 0.127, a 2 {Mz) = 0.0337. 


(28) 


The Lagrangian (|^) has terms in a form which do not allow a direct ap¬ 
plication of the non-relativistic impulse approximation for further calculation 
of the Oufdld reaction matrix element. One should rearrange the right hand 
side of eq. (^Tj) in the form of a product of two colour-singlet quark currents 
and the leptonic current. It can be accomplished by a Fierz rearrangement 
procedure and subsequent extraction of color-singlets from the product of two 
colour-triplet and colour-antitriplet quark helds. The hnal result is 


/^ALe = 2 
^eff 


(x) 



iVg + Vx){dpsJpS — -^Jt' dp fiu) + 


(29) 
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+ iVxe + v'g- Vx!)JpsJps + VNJvAJvAg] (e(l + 75)e‘'), 

The last term corresponds to the heavy Majorana neutrino contribution de¬ 
scribed by the diagram in hg. 2(a) which we also consider in the present paper 
as discussed at the end of section 2. 

The lepton number violating parameters are dehned as follows 


% 


Vx 


Vxe 


2'iias mp 

9 G%m\ rrin 


1 + 


m 


dR 


m 




\ '2 4 

7r0^2 TTlp 

6 G%m‘^ ^ rriy. 

f drt 


= 27ra'; 


\'2 

^111 

Gf< 

^ dh 


mi 

dR 


Ud) + eUu) 

4 4 


m 


dR 


, m. 




m. 




E' 

i=l 


Li\ 


m 


Xi 


Vg = 




Airas mp (mj 


9 


F^dR ^9 


m- 


7ra2 A?j 


111 


m: 


2 4 


3 Glmj 

dR 


m. 


E [^Riid)^Li{e) + 




i=l ^Xi 


Vn = 


+ eLi{u)eRi{d) 
mp 


m. 


dL 


m, 


+ (^Li{u)eLi{e) 




f^dR^ 

to 

_1 

\muRj 



< >' 


(30) 

(31) 

(32) 

(33) 

(34) 


(35) 


where < mAr > is the effective heavy Majorana neutrino mass (for dehnition 
see 0). 

Colour-singlet hadronic currents have the form 


Jps 

= h“(l-F75)da, 

(36) 

TflU 

U rp 

= h"(T^''(l-F75)d„, 

(37) 

Jav 

= h“7'"(l - 75)da, 

(38) 


where = |[7^,7‘^]- 

Now we can consider nucleon matrix elements of these quark currents which 
provide us with certain information about quark states in the nucleon. Further 
we will use in Sect. 5 the proper nuclear wave functions to describe nucleon 
states in the nucleus. Applying this standard two-step-procedure we will obtain 
the reaction matrix element for the decay. 
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4 From quark to nuclear level 


Adopting the above mentioned two-step-procedure let us write down the Oi//3/3 
decay matrix element TZoupfj corresponding to the effective Lagrangian eq. 

Using the general formula eq. (|86[) one can write 




01//3/3 


^mp^C'o^^{e(l + 75)e‘^} x 


(39) 


hlg + Vx) {p\^q\l) + {vxe + Vg- V^f) 

The heavy Majorana neutrino exchange contribution has been included in the 
last term of this equation. It corresponds to the last term of eq. 

We have introduced transition operators fli (see eq. (^) in Appendix) 
which are useful for separating the particle physics part of the calculation 
from the nuclear physics one. The transition operators contain information 
about the underlying interactions at the quark level (p^ and quark states 
inside the nucleon. They are independent of the initial |/ > and the final < F\ 
nuclear states. To calculate the nuclear matrix elements in eq. (P^D one should 
use nuclear model wave functions. 

Nuclear matrix elements of the transition operators Qq, Q? and Qn describe 
transitions induced by quark colour singlet currents (p^ ) -([3q ) in the first, sec¬ 
ond and third terms of eq. (^). Diagrams in fig. 2 with the intermediate 
states {W — N — W}, {u{d) — x, 9 — u{d)} and {u{d) — X)9~ d{u)] e — y — e, g} 
contribute to operators Qn, and flj, respectively. 

The relevant formula for calculating the transition operators in the non- 
relativistic impulse approximation (NRIA) (p^ require the nucleon matrix 
elements of these currents. 

Now we turn to the derivation of nucleon matrix elements of the colour 
singlet quark currents (P^ ) -(|35| ) using results of ref. |P3[ 
elements are, 


The relevant matrix 


< P{p)\ud\N{p') > 
< P{p)\ux5d\N{p) > 
<P{p)\ua^^’^{l + X5)d\N{p') > 
< P{p)\uY{l - X5)d\N{p) > 


Ff (g2).iV(p)r+iV(pO, (40) 

FP{q^)-N{p)x,T+N{p’), (41) 

N{p) r+iV(p'), (42) 

N{p)Y [Fv{q^) - FA{q^)l^ T+N{p%‘i) 


where q = p — p' and the tensor structure is defined as 
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Table 1: Nucleon form factor normalizations (at = 0) as calculated in 


Set 

77 ( 3 ) 

_pF) 

X p 

T{3) 

7^(3) 

7^(3) 

-^3 

A) Bag Model 

0.48 

4.41 

1.38 

-3.30 

-0.62 

B) Non-relativistic 
quark model 

0.62 

4.65 

1.45 

-1.48 

-0.66 


jrp{3} rp{3) 

^ y(3) ^ ^ (^/.P _ ^-P P) . ( 44 ) 

mp rrip 

For all form factors Fv,A{q^), we take following ref. [Q a 

dipole form 


Fv.A(q") ^ FiW) ^ Tf'if) _ d - T_ 

fv.A F^^Uo) rf’(o) V ^A 


(45) 


with rriA = 0.85GeV and fv ~ 1, /a ~ 1.261. Form factor normalization 
values Fg%{0), t/^^( 0) were calculated in ref. and are given in Table 1. 
Using formulas we may derive the non-relativistic limit mp S> |p| 


for nucleon matrix elements of the three combinations of the quark currents 
in (^). Keeping all terms up to order in the non-relativistic expansion we 
hnd from eq. (|8^) the relevant transition operators Qg, for two outgoing 

electrons in the S-wave state, 


G,- = 


T^p r ( 0 ) 


m^ 


^aP^lp^N + OiA^^GT,N + -I- a^A^GT' 


nj = = 0 ), 

^F,N — ^GT,nY 

where partial transition operators are 


^GT,N = Y. ( — I Fn{xa), 

i^j \ ^'^3 / 

^F,N = 

VijJ 


Ra 


+ (46) 

(47) 

(48) 


(49) 

(50) 
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^GT' = ■ Gj ( — ] Fi{xA), 

i^j \ / 

^F' = f—^ ^ 4 ( 0 ;^), 

T^j VijJ 

Qt’ = Y. ■ ^ii)K' • ^ij) - ■ (^j} f— ) Mxa). 

i^j \ / 


(51) 

(52) 

(53) 


Here, Rq is the nuclear radius, introduced to make the matrix elements di¬ 
mensionless (compensating factors have been absorbed into the phase space 
integrals [Q). The following notations are used 


Fp-= ( r"i - r"j), rij = \rij\, r^^ = rij/rij, XA = mArij. 


Here is the coordinate of the ”ith” nucleon. The above matrix elements 
have been written in the closure approximation, which is well satisfied due to 
the large masses of the intermediate particles. 

The nucleon structure coefficients in (|46D are given by 


a 


( 0 ) 

V 


a 


( 1 ) 

V 




fA j 


2 




jdS) I ' 

2 TP 4 j,P)\ 
3 ^ 3 j 




(54) 

(55) 

(56) 


(A'p 



1 2 \ 

3 3 j 



(57) 


Here F^% ^ F^%(0), T® ^ T®(0). 

Three different structure functions Fi appear in eqs. (iD-(0 (we use 
notations of ref. They are given by the integrals over the momentum q 

transferred between two nucleons (see eq. (|5^)) 


Fn{xa) 




1 

{2'kY {m\ + q^)^ 




(58) 
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F4{xa) 

F5{xa) 


4:7im\rij 


27im\rij 


(27r)3 (m^ + q2)4 

c;3q q2-l(q.r^^.)2 _ 
(27r)3 {rriA + 


(59) 

(60) 


These functions are analogous to the ’’neutrino potentials” for the case of 
the light Majorana neutrino exchange. Nuclear energy denominators typical 
for such momentum integrals can be neglected safely from eqs. (^^-(|60|) in 
the case of heavy intermediate particles, such as SUSY particles and the heavy 
Majorana neutrino N, with masses much larger than the characteristic nuclear 
energy and are therefore suppressed in (|^ -(|60|). The analytic solutions of the 
integrals in eqs. (|58D- (|60D are 


Fn{x) = ^(3 + 3a; + x‘^)e F 4 {x) = ^(3 + 3a; - a;^)e 


48 ' 


(61) 


F.(x) = 


XI = 


At this stage we point out that our formulas for the coefficients 
of the transition operators (E§-(0) disagree with the corresponding formulas 


derived in ref. |^. Particularly, in our case a^A — 0 while in ref. 
this coefficient is positive. This sign difference has an important consequence 
if one considers simultaneously the supersymmetric and the heavy Majorana 
neutrino contributions. It will be seen in the next section that in both the 
neutrino exchange and in the .^pSUSY mechanisms the dominant contributions 
correspond to the ^Igt,n transition operator, initiating Gamow-Teller O’*" —> 
O’*" nuclear transitions. Comparing eq. (|4^ with eq. (^8|) one can see that 
our formulas correspond to a constructive interference between the dominant 
flGT,N terms of these two mechanisms while formulas from ref. correspond 
to a destructive one. In the latter case both contributions can cancel each other 
and by a proper choice of in (^ ) -(|3^ ) the matrix element of Oi//3/3-decay 
can be set to zero at any values of particle masses involved in the formulas. As 
a result the Oi//3/3-decay half-life limit would neither constrain supersymmetric 
particle masses nor that for Majorana neutrinos. Our formulas ®), 13) and 
(|54D -(p7D always lead to certain constraints on these masses. We return to a 
detailed discussion of this point in section 6. 
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5 Nuclear matrix elements in the pn-QRPA 
approach 

While the formalism of the SUSY QuPP decay outlined in the previous sections 
is independent of the nuclear structure model used to generate the wave func¬ 
tions, numerical values of matrix elements are model dependent. We therefore 
separate this part of our work from the rest of the formalism in order to split 
up nuclear and particle physics uncertainties in a most distinctive way. 


5.1 Basic summary of the nuclear model and its param¬ 
eters 

Nuclear matrix elements have been calculated within pn-QRPA (proton - neu¬ 
tron Quasiparticle Random Phase Approximation). pn-QRPA, the RPA for 


charge-changing transitions has been developed by Halbleib and Sorenson 


and during recent years widely been applied to double beta decay calculations 


4^ - |]5^ . In the present work we follow essentially the description of |^, 


52 


Further we just briefly summarize the numerically most important features. 
We account for typically two major oscillator shells, for example Shu and 
Ahu in case of ^®Ge. Single particle energies are calculated from a Coulomb- 
corrected Woods-Saxon potential with Bohr-Mottelson parameters [^, except 


that the strength of the spin-orbit force has been reduced according to ^2 
by 9 % for a better agreement with experimental data. For the pairing 
(in Bardeen-Cooper-Schriefer (BCS) approximation) and RPA calculation we 
used consistently the Paris nucleon-nucleon potential . The strengths of the 
pairing interactions in the BCS calculation have been adjusted to reproduce 
the experimentally observed pairing gaps [^. In the RPA calculation, for the 
strength of the particle-hole interaction a simple mass number dependence was 
used, gph = 1-1-0.002A |^. For the strength of the particle-particle interaction 
numerical values of matrix elements are given for values of Qpp, which have been 


htted to experimentally measured f3^/EC decay strengths 


As mentioned in the previous sections, in the SUSY Oz//?/? decay the vir¬ 
tually exchanged particles are assumed to be heavy. Thus, the corresponding 
operators will be short-ranged. For this reason the short-range part of the 
nuclear wave functions has to be treated carefully. In addition to the hnite 
nucleon size effects, which are taken into account by the nucleon form factors 
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in momentum space (see eq. (^51) ), the nucleon-nucleon repulsion at short dis¬ 
tances must be accounted for. These short-range correlations are treated by 
multiplying the two particle wave functions by the correlation function |5S[ 


l-/(r) = l-e-“"'(l-6 r2). (62) 

The two parameters a and b can be related to each other, so that effectively 
there is only one free parameter, the so-called correlation length dehned by 

poo f) 748 

lc = - ds([l +/(r)]2 - 1) ~ , (63) 

Jo y/a 

with the standard value of Ic = 0.7 fm. 


5.2 Numerical results and discussion 


In the SUSY nuclear transition operators (|4q) - (|47|) there appear 5 basic struc¬ 
tures (|^ ) -(|53|). These, combined with the nuclear wave functions give rise to 
the following 5 nuclear structure matrix elements 


M.gt,n — {F\Qgt,n\I) — ■ cr^ f—^ Fn{xa)\I) (64) 

i^j VijJ 

Mf.n = (f = {F| Ehvy (^) FAxa)\I} (65) 

i^j 

Mgt' = {FIQgt'II) = (F| ■ a,, F.^xa)]!) (66) 

i^j \Ui/ 

Mp, = {F\ap\I) = (F| f:?^) Ft(xA)\I) (67) 

i^j \ / 

Mt' = {F\nT>\I) = (F| -f .^.)K''Ui) F5{xa)\I) 

i^j \Ui/ 

_ ( 68 ) 

^Note, that compared to the definitions of Fermi-type matrix elements for the light 
neutrino case Q , Q we took out the factor (/y// a)^ here. In case of the SUSY mechanism 
of Qu(3f3 decay all appearing combinations of coupling constants have been absorbed into 
the coefficients 
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We have calculated these hve basic matrix elements for the experimentally 
most interesting isotopes. The results are given in table 2. As expected M.gt,n 
and M.f,n have the largest numerical values, whereas the other 3 matrix el¬ 
ements give only minor corrections. This can be simply understood, if one 
notes that matrix elements A4 gt', -M-f', -Mt' have additional factors in 
the integrals (|5^- (|60D over the momentum q, transferred between two decay¬ 
ing neutrons. The effective nuclear cut-off for momentum transfer is q ^ Pf, 
where pf ~ lOOMeV is the momentum of the nucleon Fermi motion inside the 
nucleus. This then results in the relative suppression of the nuclear matrix 
elements of the last three operators in eqs. (|66|)-(|68D by factors of {pF/mpY — 
(few) %. 

The data of table 2, however, shows also another interesting fact: Corre¬ 
sponding matrix elements for different isotopes are relatively similar to each 
other, with a spread of about a factor of 2 from the mean only. As is known, 
much larger differences are usually found in calculations of 2z//3/3 decay matrix 
elements as one goes from one isotope to another. The similarity of the SUSY 
decay matrix elements is, on the other hand, not really a surprise. Recall 
that due to the large masses of the intermediate particles the transition oper¬ 
ators are short-ranged. Therefore only the part of the wave functions at short 
distances contribute appreciably to the matrix elements and nuclear structure 
effects, which are so important in 2i>/3/3 decay, play a less dominant role in 
SUSY 0i//3/3 decay. (See also the discussion in the next section.) 

According to the eqs. (@), (^3) ^ basic matrix elements (|64|)-(|6^) can 

then be combined to dehne the following nuclear matrix elements describing 
contributions to Oz//3/3 decay which correspond to the .^pSUSY mechanism 


= {F\ng\I)= (69) 

= - \cXy^^A4FN T cYjYJ iA-GT N T Ci\Y J'Af' + cYa IAgT' T 

rrie ^ 

Mj = = Wl,-|vT.=o (70) 

and the conventional mass mechanism due to heavy Majorana neutrino ex¬ 
change 

A^at = (F|ffAr|/) = M.F,N — M.GT,N^i (71) 

Table 3 shows Aiq and Ai j for the isotopes considered for the two sets of input 
parameters for the nucleon structure constants as calculated from table 1. 
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Again, it is noteworthy that both sets of parameters lead to rather similar 
resnlts. This is dne to a partial cancellation of the differences in the inpnt 
valnes for the at least in the case of Mg. 

Table 3 moreover shows that the M^s are mnch smaller than the Mg's. 
This difference can be simply understood because only scalar and pseudoscalar 
currents Js,p contribute to Mj, but no tensor current Jp- Since in this case 

(see ([70|)) a® ~ = 0, there is no contribution from the large Mgt,n 

to Mj and consequently Afj -C Mg. This result has important implications 
in the numerical analysis of section 6. 


5.3 Uncertainties of the nuclear structure matrix ele¬ 
ments 


We have investigated the dependence of the numerical values of the matrix 
elements on our choice of nuclear model parameters. Besides the numeri¬ 
cal uncertainties discussed here, there will also be deviations from the ’’true” 
matrix element due to model approximations. These, however, can not be 
quantihed exactly. Some conhdence in the model might be derived from the 
fact that the 2n/?/9 decay half-life of '^®Ge been predicted correctly 

within a factor of 2 (the 2z//d/? decay matrix element within ^/2). 

In the following the three most important parameters are discussed in some 
detail. Other parameters, like for example Qph or the use of another nucleon- 
nucleon potential (Reid soft-core potential |^) instead of the Paris force have 
been found to have negligible effects. 

In calculations of 2p(3(3 decay 


T7[| , it has been shown that the 


2v(dj3 decay matrix elements calculated within pn-QRPA do strongly depend 
on the strength of the particle-particle interaction Ppp. We have therefore cal¬ 
culated SUSY decay matrix elements as a function of Ppp for all isotopes 
considered in the present work. A typical example is shown in £g. 3. Com¬ 
pared to the results of 2n/9/9 decay calculations only a very weak dependence 
of the matrix elements on Ppp is found. A variation of Ppp within ±2 standard 
deviations from its htted value changes the SUSY matrix elements only 
by about 10 %. 


^Although Mq and are negative numerically, positive values are given in table 3. 
This was done to adjust the signs of matrix elements to the notation used in our previous 
publications on OuM decay in left-right-symmetric models Q , . 
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The explanation of this weak sensitivity of matrix elements on gpp is similar 
to the one discnssed in 0z//3/3 decay calcnlations for light nentrino exchange 
[Q. Dne to the selection rules for the Gamow-Teller operator, in 2^(3(3 decay 
the intermediate nuclear states are l’''-states, exclusively. On the other hand, 
also higher multipoles contribute to the On/3/3 matrix elements in eqs. (|6^ - 
(|68D similarly to the above mentioned case of the light neutrino exchange 
mechanism. A typical example of a multipole decomposition for M.gt,n and 
M.f,n is shown in fig. 4. Quite large contributions up to high multipoles are 
found. Since the particle-particle force mainly affects the l’''-states |^, SUSY 
On/3/3 decay matrix elements are not very sensitive to the actual choice of Qpp. 

Figure 5 shows the nuclear matrix elements for the example of ^^°Te as a 
function of the momentum cutoff factor m^. It is seen that M.gt,N) M.f,n 
and M-t' are rather stable against variations of tua, while M.gt' and M.f' are 
more sensitive. Since is dominated by A4gt,n, it will not be very sensitive 
to variations of uia as well. A4 j, on the other hand, has to be expected to 
depend more strongly on m^. 

Figure 6 then plots matrix elements as a function of the correlation length 
Ic for the example of ^®°Nd. As can be read off, a 20 % change in Ic typically 
changes calculated M.gt,n and M.f,n by about (20 — 30) %. Similar to the 
case when the cutoff factor is varied, A4gt' and A4f’ are found to be more 
sensitive to changes in Ic- 

Since changes in all three parameters can lead to both larger and smaller 
matrix elements, the influence of a simultaneous change of these can not be 
estimated by simply adding up the individual errors. To estimate the total 
uncertainty, we have therefore calculated Aiq and Ad^ at 125 grid points in this 
3-dimensional parameter space. (Parameter variations within ±20 % of their 
standard values for Ic and and within 2 standard deviations for gpp.) From 
this set we then calculated the ’’mean” matrix elements, the corresponding 
standard deviations and extreme values. 

In case oi A4q these mean matrix elements are always very similar to those 
calculated for the standard values of input parameters (typical differences up 
to a few %). The standard deviations for JHq and its extreme values deviate 
typically less than 20 % and 50 %, respectively, from the mean. 

The situation is different in case of A4 j. As discussed above, M.gt' shows 
a larger sensitivity to model parameter variations than do other matrix ele¬ 
ments. Moreover, in case of Ad ^ there occurs a destructive interference be¬ 
tween contributions from A4f,n, Mt' and Mgt'- In certain regions of the 
model parameter space M.gt' tends to cancel the contributions from the other 
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two matrix elements. The conclusion therefore unfortunately must be that J\4 j 
has to be considered to carry a large uncertainty: At extreme combinations of 
parameters Af j is smaller than (larger than) in the standard calculation by 
up to a factor of 5 (2). 

To summarize this discussion, it can be stated that while A4q can be reli¬ 
ably calculated, the numerical value of A4j must be considered to be rather 
uncertain. Fortunately we will see in the next section that only the Aig matrix 
element is relevant for the extraction of the .^pMSSM parameters from 0z//3/3 
experimental data on decay. 

6 Constraints on the ]/ip MSSM parameter 
space 

In the following we will analyse constraints on the MSSM parameter 
space using the experimental half-life limit of ^®Ge, recently measured by the 
Heidelberg-Moscow collaboration (for constraints derived from experiments 
on other isotopes, see below), 

^(^^Ge, 0+ ^ 0+) ) 5.6 x lO'^^years 90% c.l. (72) 

The theoretical expression for the inverse half-life (see (P^)) corre¬ 

sponding to the reaction matrix element TZou/sy in eq. (|39|) can be written 
as, 

[T°/f^(0+ ^ 0+)]“^ = Goi[{r]g+ri^)Mq + {r]^e+v'g-V^f)Mf + r]NMNY- (73) 

In the numerical analysis we will employ for dehniteness the matrix elements 
for Set A of the coefficients (see table 1). Matrix elements of Set B would 
only yield slightly more stringent bounds, but not change any of the arguments 
presented below. 

Neglecting for the moment the contribution from Majorana neutrino ex¬ 
change (the last term in eq. ([73|)) we will hrst analyse the supersymmetric 
contribution alone. 

Eqs. ([T2|) , (|T^) dehne an excluded area within a 5-dimensional MSSM 
parameter space 

{tan/?,/i,mg,mo, (74) 
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These parameters, as explained in section 2, completely define the snpersym- 
metric particle spectrnm and their interactions initiating Oz//3/3. However, 
there are too many free parameters to be defined from only one constraint 
(eq. Fortnnately, the analysis gets considerably simplified realizing that 

the contribntions from the five different lepton nnmber violating parameters 


rjg, rjy., riy.g, rj?, rj'- enter eq. (|73D with very different magnitndes. It can be seen 
from eqs. (pOf) - (|^) that 


Is, Vj » Vx< Vxe< Ixi ‘f 


(75) 


with the valnes of the gange conpling constants 02 and as and for any 


field composition of the nentralino states Xi (see eq. ([T8|) ). The last mass 


ineqnality in (|75|) is well satisfied even for a glnino mass rrig as large as 500 
GeV if m^.^20 GeV. The latter is guaranteed by the present experimental 
lower bound on the lightest nentralino x = Xi from the LEP experiments 
m^^20GeV. Later on we consider the glnino dominance more carefully 
and will see that it holds for rUg up to 1 TeV. Another interesting fact can 
be derived from the nuclear matrix elements M -of table 3. It can be seen 
that the value of Mg is much larger than the one of Mj^. Therefore, taking 
into account glnino dominance eq. (^), it follows that the combination rjgMg 
absolutely dominates in the half-life formula eq. (^). The advantage of this 
fact is twofold. First, only three MSSM parameters mq,mg and are 
involved in the numerical analysis. Second, as stated at the end of section 
5, the nuclear model calculations for the matrix element Mq are much more 
reliable than those for M j. Thus, we are lucky to have a theoretically well 
controlled nuclear structure dependence in the dominant term of eq. (ffSl). 
Having this in mind we proceed with the numerical analysis. 

Two limiting cases are interesting to discuss: a) = rriu^ = rrig and b) 
rriu^ S> = rriq (the case 3> rriu^ is equivalent). The following bounds 
are obtained 


Case a) < 3.9 x 


10 -^ 


Trig s 


Al/2) 

(76) 

lOOGefr^ 

' ^lOOGeV' 

j 

rUg N 

2/ m-g 5 

,(1/2) 

(77) 

lOOGeV' 

' ^lOOGeV' 

) 


Case h) A(]^^ < 5.6 x 10“^^ 

These two extreme cases differ just by a factor of 1 / 2 , and for all other ratios 
of limits between case a) and case b) are found. Motivated by the 
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MSSM mass formulas, eq. ([T^)-(p^, in the following discussion we will always 
assume = rriu^ which corresponds to the case a). 

A graphical representation of this bound (eq. ([75|)) is shown in £g. 7. The 
area on the backside of the surface is forbidden by the Oz//3/3 decay constraint 
(|7^ . Shown is the mass range between 10 GeV and 1 Tev and a coupling 
range between 10“^ and 2 x 10“^. It is interesting to note that employing the 
upper bounds rUg^g^lTeV, motivated by the SUSY naturalness argument, one 
can obtain from (|7^) an upper bound for the .^pYukawa coupling constant 


A;n~0.124. 


As mentioned before we have analysed, how the presence of a non-zero 
neutralino contribution to eq. (|73D , proportional to ^ affects the bounds 
derived from eq. 


(|761) . From eq. (^), one can expect that bounds derived 
from rjg would be sharpened when t]^ is switched on. The large difference in 
magnitude, however, leads to the result that only for very light neutralinos this 
effect will be of importance. This is displayed in fig. 8, where the excluded 
area in the (m^, rrig, m^)-space is shown for the example of = 10“^. For a 
neutralino heavier than a few GeV, the squark mass and gluino mass bounds 
get independent of the actual value of m^. As can be shown, this is true for 
any neutralino composition. Since such a light neutralino mass eigenstate is 

m 


which tell ns that > 20 


already excluded by the experiments at LEP 
GeV, the conclusion is that the limits of £g. 7 are practically independent of 
the neutralino mass. 

Now another important question arises. How much will the obtained limits 
on superparticle masses and the .^pYrkawa coupling be affected by uncertainties 
of the nuclear structure calculation? From the definition of the lepton number 
violating parameters rj eqs. o-ra and the Oz//3/3 decay constraint eqs. (|72[)- 
the following conclusion can be made. Any shift of the nuclear matrix 


element Aiq by the value AAiq, associated with the theoretical uncertainties, 
would change the extracted limits on gaugino mg^^ and sfermion niq^e 
masses as follows 


Arug^^ AMq 

- r\j - 

_ ^3.X ’ 

^Although in general all four neutralino mass eigenstates contribute to the decay 
rate, it is sufficient to consider only the lightest neutralino with mass m^. 
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rsj 


(79) 


Mu 

A'lii 

Arriq 

rriq 


rsj 


1 - (1 + 


AMq 

Mq 



Ma 



(80) 


From these equations one can see that the limits on squark masses rriq and 
the ^pYukawa coupling constant deduced from 0v(3l3 decay depend only 
very weakly on the nuclear physics uncertainties. For instance, a change of 
by even a factor of 2 changes the squark mass limit by less then 20 %. 

Up to now, we have not considered the influence of a possible contribution 
from Majorana neutrinos to Ou/3/3 decay. In this paper we bound ourselves to 
the heavy Majorana neutrino (N) contribution, as discussed at the end of sect. 
2. In section 4 it was shown that contributions to Ou/S/S decay from the SUSY 
mechanism and the ordinary neutrino mass mechanism add up coherently. 
This is opposite to the results of |ffl, where a destructive interference was 


derived. The combined constraint in the (mq,mq,(mAr))-space, where (tutv) is 
the effective heavy Majorana neutrino mass (see eq. (|35D ), is shown in £g. 9, 
for the case of = 10“^. It can be seen that limits on (ruTv) remain valid 
also, if one allows contributions from supersymmetric theories to Oz//3/3 decay 
(or, vice versa decay limits on the Ijip MSSM are not sensitive to the 

actual value of the neutrino mass). This is an important result of the present 
work, since for a destructive interference between the neutrino mass and the 
SUSY mechanisms of Oz//3/3 decay, it would always be possible to find regions 
in parameter space, where no constraints from Oz//?/? decay could be derived. 
As an interesting by-product of this analysis, we note that the current limit 
(|72D on the Oz//?/? decay half-life of ^®Ge implies 


(niTv) ~ 5.1 X 10^ GeV. 


(81) 


As mentioned above, currently ^®Ge provides the most stringent limits on 
0z//d/9 decay. For completeness, we have calculated the limits on supersym¬ 
metric parameters in the form of eq. ([75|), for all isotopes considered in the 
present work. A summary of these limits, together with the references for the 
experimental half-life limits, are given in table 4. 
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7 Comparison of constraints from 0uf3f3 decay 
with other experiments 

It is interesting, to compare limits on the MSSM parameters from Oz//?/? 
decay, with those derived from other experiments. These constraints can be 
derived from low-energy processes involving virtual superparticles or from 
direct accelerator searches Since decay is sensitive only to the hrst gen¬ 

eration lepton number violating .^pcoupling, we will restrict the discussion to 


limits on Limits on other couplings might be found in the literature []I4| , 


In ref. [15] various low-energy processes have been analysed. It was con¬ 


cluded that the most restrictive limit on might be derived from charged 


current universality. The limitation follows from the fact that the existence of 
.^pYukawa coupling X[-^LiQjDk gives an extra contribution to quark semilep- 
tonic decays (e.g., in nuclear (3 decay). The effective four-fermion interaction 
induced by the Ifip MSSM contribution in £g. 10 has a (V — A) ® {V — A) 
form identical to the one derived in the standard model. Therefore its contri¬ 
bution is equivalent to a shift in the Fermi constant Gp- If one assumes that 
only one .^poperator has a sizable coupling constant, for instance 
to a violation of charged current (CC) universality II because the shift of 
Gf is different for different generations. The experimental limits impose the 
following bound at the two sigma level 


A'ni < 


0.03f- 


m 


dR \ 


lOOGeV^' 

If one allows more than one ^p operator to contribute the violation of CC 
universality is reduced and the above bound is weakened. Comparing (|^ and 
(|8^ one can see that for masses in the range of 100 GeV, the bound (^) is 
less restrictive by nearly 2 orders of magnitude than the bound ([76| ) derived 
from Oz//3/3 decay. 

Accelerator searches on supersymmetric particles are usually based on the 
assumption of R-parity conservation. In this case the LSP, which is assumed 
to be the neutralino X; is stable and only weakly interacting. Therefore it 
escapes from the detector yielding the prominent missing transverse energy 
{l^T ) signature of SUSY events. The CDF bounds on superparticle masses 
rely essentially on .^^signals and are not valid in the ljipca.se. However, in 
ref. |T^ it has been shown that limits on superparticle masses in the 1/ipCase 
can be derived using the CDF dilepton search data. The corresponding mass 
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limit is 


'^g,q > lOOGeV, (83) 

independently of if > 10“^ is assumed. For smaller values of the 

LSP has a negligible decay probability inside the detector and the signal 
has pattern as in the Rp conserving case. Then the limits obtained in ref. 
^ can be applied for such small .^pcouplings. A recent calculation shows 
that even larger masses than in eq. (|83|) might be probed at the TEVATRON 
in the near future. 

There are in the literature other proposals for searching .^pSUSY signals 
in future accelerator experiments. 

The isolated like-sign dilepton signature is proposed as a characteristic 
feature of the .^pevents at the LHC 

Recently searching for .^pSUSY events in deep inelastic ep-scattering ex¬ 
periments with the ZEUS detector at HERA has been proposed [^]. Two 


sets of signals could be identihed with these events. They are the .^presonant 
squark production followed by the .^pcascade decay of the neutralino and the 
MSSM Rp conserving production of selectron and squark with their subsequent 
.^pdecay to ordinary matter. It was advocated that searching for these signals 
provide HERA with a promising discovery potential for .^pSUSY. 

A comparison of the bounds, which according to might be reached with 
one year of HERA data, and the other limits discussed above with the limits 
we have obtained from the decay constraint (^) is shown in £g. 10 in 

the — rriq plane. In the case of QuPP decay, limits for 2 different values of 
the gluino mass are shown. It can be seen, that even for a gluino mass as large 
as 1 TeV which is marginal from the point of view of SUSY naturalness, the 
present double beta decay half-life limit yields the most restrictive bounds. | 

To summarize this discussion, Oi//3/3 decay allows to stringently restrict Rp 
violating supersymmetric theories. We have shown that these limits are more 
stringent than those from other low-energy processes as well as those which 
can be derived from the HERA experiments (see fig. 10). 


8 Conclusions 

We have investigated contributions to 0z//3/3 decay from supersymmetric theo¬ 
ries with explicit R-parity violation. The complete set of diagrams describing 

'^This result, however, does not touch upon other couplings than which are 
unaccesible for 0vf3(3, but can be probed in other processes. 
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quark-lepton interactions at short-distances has been considered. On this ba¬ 
sis we have obtained the relevant low-energy effective Lagrangian in terms of 
nucleon and lepton currents. Then we have derived the transition operators 
describing O’*" —O’*" nuclear transitions induced by the .^pMSSM interactions. 

It has been found that contributions from the .^pSUSY mechanism of QuPP 
decay add coherently to the well-known neutrino mass mechanism. As a result, 
limits on the Majorana neutrino mass derived from QuPP decay remain valid 
also if the supersymmetric contributions are taken into account. 

We have calculated nuclear matrix elements of these transition operators 
within a realistic nuclear structure model for experimentally interesting iso¬ 
topes. 

Special attention has been paid to the theoretical uncertainties of the nu¬ 
clear matrix elements and to the question how these uncertainties affects limits 
on the MSSM model parameters we extracted from Qu(3(3. We were able 
to conclude that these limits do only very weakly depend on nuclear physics 
uncertainties. 

Using existing experimental lower bounds on the 0z//5/5 decay half-life T 1 / 2 , 
we then analysed constraints on the .^pMSSM parameters from 0z//5/5 decay. 
The most restrictive limits are found from the current ^®Ge QuPP decay ex¬ 
periment by the Heidelberg-Moscow collaboration. We conclude that QuPP 
decay imposes very restrictive bounds on the lepton number violating sector 
of .^pSUSY models. We have presented these bounds as a 3-dimensional exclu¬ 
sion plot in the space of the .^pYukawa coupling constant squark rriq and 
gluino rriq masses (fig. 7) as well as a 2-dimensional one in the X'^-rriq plane 
(fig. 10) which also shows bounds from other low and high energy processes. 

We infer that the Qvj3j3 bounds are able to compete with or are even more 
stringent than those derived from current and near future accelerator exper¬ 
iments. Particularly, they exclude the domain which is accessible for the ex¬ 
periments with the ZEUS detector at HERA. 
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9 Appendix 


Here we define some notations used in the main text. A complete set of the 
relevant formulas can be found in ref. P]. 

The 0z//3/3 reaction matrix element TZou/b/b is dehned as follows 


(FIS' - 1|/) = 27ri6{ei + 62 + Ef - Ei)Tloy 0 fi (84) 

where Fj, Fj, ei ^2 = energies of the initial and hnal 

nuclear states as well as the energies of two outgoing electrons with the 3- 
momenta p). 

For the Lagrangian represented in the form 

F ={e(l+ 75)e‘'}^PiJiJi, (85) 

i 

the corresponding reaction matrix element can be written as 

{e(l + 75 )e‘^}^pi(F|f]i|/) (86) 


The only approximation made here is that nuclear matrix elements are inde¬ 
pendent of the hnal-state electron energies 61^2 and momenta pi^ 2 - This ap¬ 
proximation is well satished in the Oi//3/5 decay, because ei, 2 ~To ~ (1 — 3)MeV 
(To is the energy release) while the typical energy scale of the nuclear matrix 
elements is given by the nucleon Fermi momentum, pp ~ lOOMeV. Corrections 
to eq. (| 8 ^ are smaller than 1 %. 

We introduce transition operators ffj as 

(F|a|/) = 7 ^ jd\Ff(<f) (87) 


The reaction matrix element (j^) does not depend on the value of the numer¬ 
ical coefficient which is introduced to bring the normalization of hi in a 
coincidence with the literature 


Cc)y = dvr 


mp Rq _2 
- -n^A 

me fi 


( 88 ) 


The formula (|87D is written for the case of heavy intermediate particles, in the 
closure approximation and for the outgoing electrons in S-wave states. 







It is a common practice to use the non-relativistic impulse approximation 
(NRIA) for nuclear matrix element calculations. In this approximation the 
transition operator fl,, describing NaNf, —>■ PaPb transitions of two initial 
neutrons into two hnal protons induced by the current Jj, takes the form 




(89) 

^ / a^b'' 

X < PaiPa)\Ji\Na{p'a) >nr< A^b(pfe) | | A^b(p'fc) >nr ■ 


Here q = Pa ~ Pa ~ Pfe ~ Pfc momentum transfer, Ri(q^) are the nu¬ 

cleon form factors and the subscript | >nr means non-relativistic limit for the 
corresponding nucleon matrix element. The summation over all pairs NaNb of 
initial neutrons is implied. Tab is the separation between these two neutrons. 

The Oz//3/3 observables can be calculated on the basis of the reaction matrix 
element TZoupid- 

The differential width of Oz//3/3 is 


dr — 27rd(ei -|- 62 + Ef — Ei)\TZQ^f}i^\‘^d^lid^l2^ (90) 


the phase space factors are 


dflj = 


d^Pi 

(27r)3' 


The 0z//3/3 half-life formula can be written as 


^ 1/2 




(91) 


(92) 


where Gqi is the leptonic phase space integral, calculated according to the 
prescription of 0 and has the form 


Gni = 




[meRo)‘^ln{2) 


d^id^2^oi • 


(93) 


Here, ao^y = (Gi7’Cos(0c)/A//y)^"ig/(647r^h) involving only physical constants, 
Ro is the nuclear radius. The kinematical factor 601 accounts for the Coulomb 
distortion of the electron waves and can be taken from ref. 0 . 
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Figure Captions 


Fig. la Basic diagram for neutrinoless double beta decay. 

Fig. lb Basic diagram for neutrinoless double beta decay at the quark level. 

Fig.2a Feynman graphs for the conventional mechanism of decay by mas¬ 
sive Majorana neutrino exchange. 

Fig.2b Feynman graphs for the supersymmetric contributions to decay. 


Fig.2c New ’’non-diagonal” Feynman graphs for supersymmetric contributions 
to 0z//3/3 decay [^ 


Fig.3a SUSY 0z//3/3 decay matrix elements as a function of the particle-particle 
strength Qpp, for the example of ^®Ce. Shown are Mgt,n (full line) and 
(dashed line). 

Fig.3b As hg. 3.a, but for M.gt' (full line) and M.pi (dashed line) and M.t' 
(dash-dotted line). 


Fig.4 Multipole decomposition of the matrix elements A4gt,n and M.f,n for 
the example of ^®Ce. 
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Fig.Sa 
Fig. 5b 
Fig.Ga 
Fig. 6b 
Fig. 7 


Fig.8 
Fig.9 
Fig.lO 
Fig.11 


Matrix elements M.gt,n (full line) and M.f,n (dashed line) for ^^°Te as 
a function of the momentum cutoff factor ruA [MeV]. 

As £g. 5.a, but for A4gt' (full line) and A4f' (dashed line) and M.t' 
(dash-dotted line). 

Matrix elements M.gt,n (full line) and M.f,n (dashed line) for ^®°Nd as 
a function of the correlation length Ig [fm]. 

As fig. 6.a, but for A4gt' (full line) and A4f' (dashed line) and M.t' 
(dash-dotted line). 

Constraints from Oi//3/3 decay on the squark rriq and gluino trig masses 
and the R-parity violating Yukawa coupling constant Values on the 
backside of the surface are forbidden by non-observation of Of(3(3 decay. 

~ '^UL ~ '^9 assumed and the matrix element Aig for Set A of 
the coefficients has been used. 


Dependence of the squark and gluino mass limits derived from 0Ff3f3 
decay on the actual value of the neutralino mass for = 10“^. 

Combined limits on supersymmetric particle masses rriq^g and the effec¬ 
tive heavy Majorana neutrino mass {rriN) for A(]^^ = 10“^. 

Feynman graph for neutron decay in Rp-violating supersymmetric theo¬ 
ries, see text. 


Comparison of limits on R-parity violating supersymmetric theories from 
different experiments in the (mg-A(^]^) plane. The dashed line is the limit 


from charged-current universality according to S. The vertical line is 


the lower limit on squark masses in Rp-violating supersymmetric theories 
from the TEVATRON, according to |T^. The thick full line is the region 
which might be explored by HERA with about one year of data . The 
2 full lines to the right are the limits obtained from non-observation of 
the 0^(313 decay of ^®Ge for gluino masses of (from left to right) rrig = 1 
TeV, 100 GeV, respectively. The parameter regions to the upper left of 
the lines are forbidden by the different experiments. It is seen that even 
for a gluino mass as large as 1 TeV, Of(3(3 decay gives the most stringent 
limits. 
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Table 2: Nuclear matrix elements for SUSY 0u(3(3 decay for the experimentally 
most interesting isotopes calculated within pn-QRPA. xlO* implies that 
the matrix element should be divided by 10* to get the correct numerical value. 
These factors are introduced here, just to show the differences in magnitude 
of the different matrix elements. 


Ay 

M.gt,n xlO^ 

A4f,n xlO^ 

■M.QT' xlO^ 

Mf' xlO^ 

X 

o 

CO 

^‘^Ge 

1.13 

-4.07 

-7.70 

3.06 

-3.09 

82Se 

1.02 

-3.60 

-7.13 

2.76 

-2.76 

looMo 

1.29 

-4.89 

-7.88 

3.52 

-4.93 

ii^Cd 

0.75 

-2.71 

-5.49 

2.19 

-2.65 

i 28 xe 

1.19 

-4.19 

-7.95 

3.11 

-3.85 

laoTe 

1.05 

-3.69 

-6.97 

2.73 

-3.53 

i36Xe 

0.58 

-2.03 

-3.81 

1.49 

-1.81 

150Nd 

1.65 

-5.91 

-10.4 

4.32 

-7.51 
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Table 3: Nuclear matrix elements for SUSY 0i/(3(3 decay. Shown are Aiq and 
A4 j for the two sets of input values of coefficients for the of table 1. 


Ay 

A) Mq 

B) Mq 

A) Mj 

B) Mj 

^‘^Ge 

283 

304 

13.2 

20.7 

82 Se 

253 

272 

11.3 

17.9 

looMo 

328 

356 

23.6 

33.5 

ii^Cd 

190 

205 

11.0 

16.2 

i28Te 

298 

323 

16.7 

24.8 

i30xe 

262 

284 

15.4 

22.5 

i36Xe 

143 

155 

7.91 

11.8 

150Nd 

416 

456 

34.4 

47.0 
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Table 4: Comparison of limits on supersymmetric parameters derived from dif¬ 
ferent (5(3 decay experiments (see quoted references for half-life limits), written 
in the form of eq. ([76|): 


'm<4T5SfcF) ( 


lOOGeVj 


Currently ™Ge provides the most stringent limits. 
[Ref.] ^ 68 % c.l. only.) 


[Ref.] “ 90 % c.h, whereas 


Ay 

e 

Ref. 

Ay 

e 

Ref. 

^^Ge 

3.9 X 10-^ 


128yg 

4.9 X 10-^ 


60 

b 

82 Se 

1.1 X 10-8 


61 

b 

180yg 

1.1 X 10-8 


32 

a 

looMo 

7.5 X 10-^ 


63 

b 

886 Xe 

6.8 X 10-^ 


34 

a 

ii^Cd 

2.5 X 10-8 


35 

a 

150Nd 

9.7 X 10-^ 


36 
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Figure l.a 
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Figure 6.a 
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